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Abstract 

Recently, analytical solutions of a nonlinear Fokker-Planck equation describ- 
ing anomalous diffusion with an external linear force were found using a non 
extensive tliermostatistical Ansatz. We have extended these solutions to the 
case when an homogeneous absorption process is also present. Some pecu- 
liar aspects of the interrelation between the deterministic force, the nonlinear 
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diffusion and the absorption process are discussed. 
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I. INTRODUCTION 

The ubiquity of the anomalous diffusion phenomenon in nature has attracted the interest 
of researchers from both the theoretical and experimental point of view. Anomalous diffusion 
has been found in transport of fluids in porous media and surface growth j^, in NMR 
relaxometry of liquids in porous glasses 0, in a two dimensional fluid flow [^, to name 
just a few among the large variety of physical phenomena where it is present. A related 
aspect is the case of density dependent diffusivities as found in some biological systems [Q , 
in polymers 0, and hydrogen diffusion in metals (see also 0]). 

Some recent papers have investigated a class of nonlinear generalizations diffusion and 



Fokker-Planck equations flpHlOl, as a model of correlated anomalous diffusion. Some of 



those studies were based in a nonextensive thermodynamical formalism ||Tl|. Particularly in 



Ref. |]TOl exact solutions for the nonlinear Fokker-Planck equation subject to a linear force 
have been found. Here we want to show how these solutions can be extended to the case 
where an absorption process is also present. 

We start recalling the "full" anomalous diffusion equation or "nonlinear" Fokker-Planck 



equation solved in Ref. |10 



^^[pix,t)r = -^{FixMx,t)r} + D^[pix,t)Y. (1) 

When F{x) = 0, Eq. (P can be interpreted as a diffusion equation for $(x,t) = [p{x,t)]'^, 
where the diffusivity depends on $ 

(2) 
(3) 

There are several real situations where this power-law dependence of the diffusivity is found. 
It occurs in the flow of gases through porous media {i^/fi > 2 [|12|), the flow of water 
in unsaturated soils {i^/fJ^ = 5 fl^), the simultaneous diffusion and adsorption in porous 
samples where the adsorption isotherm is of power-law type []T1| (z^//U > 1 for Freundlich type 



dt 


dx^ 


(Dim) 




/}($) 





of adsorption isotherm |T4| , p!5| ). Clearly, in those cases, the diffusivity vanishes (diverges) 
for $ = when z///i > 1 (z^//i < 1). It is worth to remember that Eq. (|T]) corresponds 
to the so called "Porous Media Equation" when fi = 1 P,p^. There are a large number 
of situations where the interest of describing anomalous diffusion plus absorption are of 
relevance. Notably, those related with diffusion of some (reactive) substance in gaseous 
phase through a porous media or a membrane, that react and can be adsorbed in sites 



inside the pore [|17 



Our interest here is to solve the same Eq. (||) but including now terms that describe 
some kind of absorption process. A general form of such an equation is 

^^[p(x,t)r = -^{F{xMx,t)r}+D^[p{x,t)r-a[p{x,t)r\ (4) 

where a plays the role of an absorption rate (and becomes the usual one for fi' = fi). The 
presence of reaction terms like the one in Eq. (^ (with a 7^ and fi' 7^ 0) is not at all 
unexpected considering the large amount of work on the problematic of diffusion-limited 
reactions. Among the diversity of systems that have been studied we only recall here the so 
called one species coagulation, that is: A+A -^0 or mA—^lA (with m > I), that have been 
associated, among others, with catalytic processes in regular, heterogeneous or disordered 



systems |jT8|. The reaction term may account, in the case /i = /x', for an irreversible first- 
order reaction of the transported substance so that the rate of removal is a C fl^ . This 
extra term also appears when a tracer undergoing radioactive decay is transported through 
a porous medium, where a is the reciprocal of the tracer's mean lifetime |jl9[, as well as in 
heat flow involving heat production at a rate which is a linear function of the temperature 
pO| . Finally, in solute transport through adsorbent samples the adsorption rate, at small 
solute concentration, is usually proportional to the concentration in solution and Eq. (^ 
applies. 

In |]TU[ it has been shown that Pq{x,t), the solution of Eq. (|I]), for a linear force F{x) 
has the form 
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where q = 1 + fi — u, P{t) depends on the width of the distribution, XMif) is the average 
of the coordinate and Zq{t) is a normahzation factor. All of them depend on the diffusion 
parameter D as well as on yu, v and the force (see |T^ for details). 

For completeness, as well as for reference, we write Eq. (|1|) without the "drift" F{x) as 
we will refer to the solutions of such an equation in the following sections 

|[p(x,tr = D^[p(x,t)r. (6) 

According to the results from Ref . [|lOl we can write that the solution p^°) (x, t) of Eq. (^ 



has the form given in Eq. (^), with 

XM{t) = Xq 



/ 2z/ \ "+'' 

Zg{t)= i—{u + ix)nDt\ 

f3{t) = rf (j{u + fi)nDt\ "^\ (7) 

where, we have used the relation /?(0)Zg(0)^'^ = tt, that shall be fulfilled if we want to have 
a (5-like initial condition. 

In the present work we intent to analyze the specific case of a linear drift, namely 
F{x) = ki — k2X. This case, where the potential is harmonic (a typical approximation), is 
the simple nontrivial one where analytic solutions can be obtained just by means of changing 
the variables to suitable ones, namely a simple extension of the well known Boltzmann 
Transformation ||21|| . 

In section |I| we start considering the simple case fi = fi'. We analize first the case of 
constant external force {k2 = 0). In this case we can firstly reduce Eq. (Q) by proposing a 
solution of the form 

p{x,t) = e^'i^ p{x,t), (8) 

that yields an equation for p{x, t) given by 

^^WM = -§^{n^)W.t)Y} (9) 



Although this reduction to a nonlinear Fokker-Planck equation hke in Eq. (|l|) looks to be 
always possible (when ^2 = 0), for the case of a linear force (when ^2 7^ 0) such a reduction 
is not possible to make in a simple way and we will need a more general treatment. 

The linear force situation, tightly related to the so called Uhlenbeck-Ornstein process 
(fci = 0; A;2 7^ 0), is treated in section |ITT|, while in section |V| we discuss the most general 
case, that is when /i ^ fi'. In the last section we make some final remarks. 

II. SOLUTION FOR A CONSTANT FORCE 

As indicated above, here we consider the case of a constant force, that is F{x) = ki. 
Equation (P) can be further reduced making the following change of variables ^ = x — kit, 
that results in 

|[p(e,tr = De"(i-^/^)*^[p(e,t)]^ (10) 

Now we change the time variable according to 

p(e,t)-P(e,^(t))^| = i(t)^ (11) 

^(t) = /* e-(i-/^)-rfr = ^^^ t>0, (12) 



/o 7 

with 7 = —a{fi — i^)//i, and obtain the following equation valid for all t > 

§;m,^)r = D^,m,z)r. (13) 

Hence, if p^°^(a;, t) is the solution of Eq. (^, the solution with F{x) = ki plus absorption 
results to be 



Pi 



ix,t) = e-'^'pf\x-kit,z{t)). (14) 



It is easy to check that this solution has the right limits for a -^ and for a > and 
H = u = 1, i.e the standard Fokker-Planck equation (plus absorption). 

The new variable z{t) plays the role of an effective time for the dispersion process. It 
exhibits markedly different behaviors depending on the ratio fi/u. The case fi/u > 1 (7 < 0) 

6 



corresponds to superdiffusive transport when there is no absorption [|T0|. When absorption 



is present this superdiffusion is enhanced since the effective time z grows exponentially as a 
function of the real time t. 

In the case fi/u < 1 (7 > 0), which leads to suhdiffusion for a = 0, the presence of 
absorption also plays a key role. The effective time z{t) converges to an asymptotic value: 
limt^oo z{t) = Zoo = 1/7- Therefore, the distribution p^°''(^, z) evolves toward an asymptotic 
curve p^°^(^,Zoo)- Let us note that z^o diverges whenever a ^ or fi/u -^ 1. In Fig. 
ID we compare, in the /i/z/ < 1 case, the time evolution of the distributions for a 7^ 
and a = 0. We also compare these distributions with the shape of the a 7^ asymptotic 
curve. For completeness, we include in Fig.^ the behavior of z on t, also illustrating its 
dependence on /i/z/. Finally in the case ^ = v {normal diffusion+ absorption ior the quantity 
<l>(x,t) = [p{x,t)]^) the change of variables becomes a simple time scaling. 

As already mentioned, Eq. (^ can be viewed as a classical diffusion equation for $(a;, t) = 
[p(x,t)]'', where the diffusivity depends on $(a;,t) through -D($) = D^'^^^^^^. Therefore, it 
becomes clear that the absorption can enhance (reduce) the diffusive transport whenever 
fi/u > 1 {fi/i^ < 1), namely: as absorption proceeds $ decreases, yielding to an increase or 
not in -D($) according to the ^/u ratio. This qualitative description seems to be in complete 
agreement with the previous quantitative results. 

III. SOLUTION FOR A LINEAR FORCE CASE 

We now consider the case of a linear force, given by F{x) = ki — k2 x (here and in what 
follows we assume that k2 > 0), whose general solution without absorption was found in 
||1CI|| . To start with, we assume that /i = fi'. With the hint of the change of variables made 
in the previous section we propose the following changes, that define the new variables 

pix,t)=e-'"'pi^,zit)) (15) 

C = xg{t) + f{t), (16) 



with w, g{t) and f{t) to be determined. In terms of these new variables the time and space 
derivatives becomes 

|4 + (.«W + /M)| + .(4 (17) 

Taking into account these results and the form proposed for the solution (given by Eq. 
(|15|)), each separate term of Eq. (|) becomes 

^^[pix,t)r = -wf,e-'"^'m,z)r + (19) 

+ e-'"^*i(t)^[p(e,;^)r, 

- ^ {{h - k2xMx,t)r} = k2[p{x,t)r - ih - hx)^[p{x,t)r (20) 

D^2[pi^^t)r = Dg\t)e—'^,m,z)r. (21) 

In this way, the equation we obtain for p{^, z) replacing into Eq. (^ with F{x) = ki — k2X 
results (after arranging terms and multiplying by e'"^*) 

z{t)-^[p{^,z)r = [wix + k2-a][p{^,z)r (22) 

-[(fci - k2x)g{t) + X git) + /(t)]^[p(e, z)Y 

+D/(t)e-(^-'^)*^[p(e,^)]^ 

In order to reduce the last equation to a one with a form similar to Eq. (|^), we need to cancel 
the first two terms on the rhs, and reduce the coefficient of the third one to a constant. To 
operate with the second term, we shall cancel it for all values of x. These conditions yield 
the following equations 
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= wjji + k2 — a (23) 

= -k2g{t) + g{t) (24) 

= hg{t) + fit) (25) 

1 = g\t)e'"^^'-''^'z{ty\ (26) 



rendering 



w = {-k2 + a)ij-^ (27) 

g{t) = Ge'-' (28) 

fit) = H- G^e^'' (29) 

zit) - ziO) = G^ {l - e'^'} ^~\ (30) 

with 7 = i—k2ifi + z/) + q;(z/ — fi)) fi^^. In the general case, the values of the constants shall 
be chosen to fulfill some particular initial condition. Here, to simplify, we choose G = 1, 
implying that we do not change the x scale at t = 0. Also, in order to make the change of 
space variables in such a way to have it centered at the potential minimum {(,' = (x — |^)) 
we adopt H = 0. Finally we choose z(0) = to preserve the time origin. With these values 
we have 

^^(.-|)e- (31) 

z(t) ^ {l - e--"} j-\ (32) 

and the solution of Eq. (^ with F{x) = ki — k2 x is 

p,ix,t) = e^'pf{i,z). (33) 

This is the final result for the present case (yU = /i'). It is trivial to check its validity in 
some limits, the most obvious one is to choose a = 0, recovering the solution of Ref. []10[ . 



With yU = z/ = 1 and a 7^ we recover the simple case of diffusion in a harmonic potential 
with absorption. Also, if we consider the case ^ = v and a 7^ 0, it is immediate to obtain 
(remember that ^ = v gives g = 1!) 



This result becomes obvious after making the change pi{x,t)^ = (f){x,t), reducing the prob- 
lem to an effective one of diffusion in a harmonic potential with absorption for (j){x,t). 
Clearly, even though the solution has a Gaussian form (times a decaying exponential term), 
the width of the Gaussian factor behaves "anomalously" as it differs from the one in the 
associated Ornstein-Uhlenbeck process P^ . 

As in the constant-force case, absorption process markedly influences the time evolution 
of dispersion. A straightforward calculation yields the dispersion of the distribution in the 
present case (i.e. fj, = fj,') 

In the super diffusive case (yu/z^ > 1) P{z{t)) becomes asymptotically exponential, namely, 
P{z{t)) (X exp (—2 ^'■^+'^-'+"'-^~'^^ tj. Replacing this result in Eq. (|35|) we obtain the long time 
behavior of the dispersion 

((x-(x))2) oce'"(^*. (36) 

Therefore, the superdijfusive transport enhanced by absorption yields an exponentially in- 
creasing dispersion even in an attractive potential. The subdijfusive case (yu/z^ < 1) presents 
two different situations. Although in both cases the dispersion decays exponentially, when 
the absorption rate is small (7 < 0, a < k2{i^ + ft) / {i^ — ft)) , absorption is the rate controlling 
process for dispersion 

{{x - (x))2) oc e-'"fe^*. (37) 

In the other case, when the absorption rate is large (7 > 0), the attractive force becomes 
rate limiting for dispersion process 

{{x - (x))') oc e-^'''\ (38) 
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In order to compare the influence of tfie absorption term on tlie solutions we have found, 
in Fig. ^ we depict the solution given in Eq. (p^), in the case ^/v > 1 (superdiffusion), for 
a = and a 7^ at different times. In Fig. ^ we compare, in the subdiffusive case, the 
solution in Eq. (|33| ) when absorption is the rate limiting process for dispersion to the case 
when the attractive force controls dispersion. In both figures the differences between the 
characteristics of the different situations are apparent. 

IV. GENERAL ABSORPTION TERM 

In this section we consider Eq. (|), in the general case fj, ^ fj,'. The following "simple" 
kinetic equation 

^^[p{t)r = -a[p{t)f, (39) 

whose solution is 

a 1 

p(t) = [l_(l_g')_t]-(W), (40) 

A* 
where g' = 1 — /i' + /i, strongly suggests to replace the exponential in the change of variables 
in Eqs. (|T^ and (|T3p, by the q'-exponentialiunction defined by Eq. (^Op [^. The ordinary 
exponential function is recovered when q' —>■ 1. If we try this possibility, together with 
the Ansatz in Eq. (^, it immediately leads to the condition /i = /i'. This result becomes 
apparent when analyzing Eq. (Q) in terms of $(x, t) = [p{x,t)]'^. The form proposed in 
Eq. (^ allows to reduce the general equation, eliminating the absorption term, only when 
absorption is proportional to $(a;,t). However, the general case with fi ^ fi' will have a 
solution whose scaling properties can be determined. 

To find such scaling behavior it is enough to consider the simplified situation without 
external force, that is 

|[p(x,t)]'^ = ^^b(^'^)]^ - ^[pi^^t)f- (41) 

We consider the following Ansatz 
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p{x,t) = ip{t)e{0, 



(42) 



with ^ = ^(t) X. Replacing this into Eq. (PTI), we obtain the functions ip{t) and iplt) as 



^(t) = [1 + (1 - g') (a//i)t] 

V>(t) = [1 + (1 - g') (a//i)t]-^(W) 



/ 
1 fl —u 



(43) 
(44) 



where q' = 1 — fi' + fi. Hence, Eq. (^) for 0(0 reduces to an ordinary differential equation 
on the variable ^ 



e^ 



jj, — u 
2/i 



«|«'' = "|^«"-^®'" 



(45) 



Once again, the previous results can be interpreted in terms of $(x,t). Writing the ab- 
sorption term as —{a^^ /^^^)$, it can be seen that in the case /i' < /i the absorption process 
is enhanced as $ decreases with time. This leads to a finite time tc = fi/a{^ — /x'), where 
$ becomes zero. On the other hand, when /i' > /i we can obtain the asymptotic dispersion, 
even though the ordinary differential equation for 0(0, Eq. (|i5|), is too complicated to be 
solved analytically. However, for the nt/i-moment of the distribution, we obtain 



(x-) 



dx x^" p{x, t) 

dx x^" ip{t) 0(V^(t) X 
V'(t)-2" [ / d^ e^ 0(0 
V'(i)-'M2„ 



dx p{x,t) 

dx ip(t) e{ij{t) x) 
d^ 0(0 



(46) 



{x'^+^) = 0, 



(47) 



yielding 



((x - {x)y) = (x^) = ^it)-^A2 ~ tV(i-«') 






(48) 



Hence, it is clear that, as in the previous case (/i' = /i), fi/u < 1 corresponds to subdiffusion 
whereas the case fi/u > 1 corresponds to superdiffusive transport. 
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V. FINAL REMARKS 

The Fokker-Planck equation that was generahzed to a nonextensive scenario has been 
further generahzed to include the possibihty of an absorption process. We have shown that 
the exact solutions of Eq. (^ (a nonlinear Fokker-Planck equation subject to linear forces) 
found in [|I^ when a = 0, can be extended for the case a 7^ and /i' = fi. However, 
in the general case fi' 7^ /i, we have been only able to obtain the scaling properties of the 
solution (whose analytical form we cannot obtain), and the asymptotic behavior of the whole 
hierarchy of moments. 

Summarizing our results for the nonlinear process of anomalous diffusion plus absorp- 
tion, as described by Eq.(|), we have found that the solution 

- in a constant force field and for fi/u > 1 shows a super diffusive behaviour that is enhanced 
when a 7^ (7 < 0), 

- also, when 7 > (/i/z^ < 1) the concentration reaches an asymptotic constant profile, 

- for a linear force and 7 < 0, we find an exponentially increasing dispersion for superdiffu- 
sion, 

- also, in the linear force case, an exponentially decreasing dispersion arises for subdiffusion 
(gamma i 0), where absorption is the rate controlling process for dispersion when absorption 
rate is small {a < k2{i' + fi) / {ly — fi)) , while the attractive force becomes the rate limiting 
dispersion process when absorption is large enough (a > k2{i' + n) / {u — jj)) ■ 

The present results gives further support to the argument llO| that a generalized ther- 



mostatistics including nonextensivity constitutes an adequate framework within which it 



is possible to unify both normal and correlated anomalous diffusion ||I0|, extended now to 
the case when an absorption process is also present. Also, as indicated in 0, this kind of 
work points out the convenience of paying more attention to the thermodynamic aspects 



of non-Fickian diffusion. Moreover, it has been suggested [jTO| that even Levy-like anoma- 
lous diffusion |Q (that can be discussed by means of linear Fokker-Planck equations with 
fractional derivatives) can be included within the present common framework of nonlinear 
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Fokker-Planck equations with fractional derivatives. 

This work also opens the possibility of analyzing reaction-diffusion systems on a fractal 
substratum, by considering nonlinear Fokker-Planck equation with other forms of reaction 
terms. This problem will be the subject of further work. 
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FIGURES 
FIG. 1. The evolution of the distribution Pq {C, z{t)), as given in Eq. (0), in the case 
/i/i/ = 2/3 < 1. It is shown at times t = 0.2, 2.5 and 5 for a = and a = 1 
(/u = l,z^ = 1.5, D = 1,/ci = 5,q = 0.5). The dashed hue corresponds to the a = 1 asymp- 
totic distribution shape. 

FIG. 2. Time dependence of the effective time z for several values of 7. It is apparent that z{t) 
saturates at a finite value as i — > cxd (zoq = I/7) whenever 7 > 0. In the case 7 = the effective 
time becomes equal to the real time: z = t. For 7 < the effective time grows exponentially with 
t 

FIG. 3. Evolution of the distribution Pq {S,, z{t)), as given in Eq. (|33|), in the case 



fi/u = 2 > 1. It is shown at times t = 0.5,0.8 and 1.0 for a = and a = 5 
{fj, = IjV = 0.5, D = l,ki = 0,k2 = 1.0, g = 1.5). The dotted lines correspond to the a = 
distribution, while the solid lines correspond to the a = 5 case. 

FIG. 4. Evolution oi pq {S,,z{t)), as given by Eq. (|3^) , in the case fi/u = 2/3 < 1. We show it 
at times t = 0.2, 1 and 5 for a = 0.5 and a = 2.0 (/u = 1, z^ = 1.5, D = l,ki = 0,k2 = 0.2, q = 0.5). 
The dashed lines correspond to the case a = 0.5 where the attractive force becomes dispersion rate 
limiting. The solid lines correspond to a = 2.0 where absorption is the rate controlling process for 
dispersion. 
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